1 



Fourth-order perturbative equations in Lagrangian perturbation 
theory for a cosmological dust fluid 

Takayuki Tatekawa^'^ *) 

^ Center for Information Initiative, University of Fukui, 3-9-1 Bunkyo, Fukui 

910-8507, Japan 

^Research Institute for Science and Engineering, Waseda University, 3-4-1 Okubo, 
CN " Shinjuku-ku, Tokyo 169-8555, Japan 

, We have derived fourth-order perturbative equations in Lagrangian perturbation the- 

■ ory for a cosmological dust fluid. These equations are derived under the supposition of 

■ Newtonian cosmology in the Friedmann-Lemaitre-Robertson- Walker Universe model. Even 
if we consider the longitudinal mode in the first-order perturbation, the transverse mode 
appears in the third-order perturbation. Furthermore, in this case, six longitudinal-mode 
equations and four transverse-mode equations appear in the fourth-order perturbation. The 

' application of the fourth-order perturbation leads to a precise prediction of the large-scale 

structure. 
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§1. Introduction 



Describing the large-scale structure in the Universe is one of the most impor- 
tant problems in modern cosmology. It has been postulated that this structure was 
^ I formed from the primordial density fluctuation due to self-gravitational instability 

[ij,[21i,|3li,|4),[g,[6li,[3,|8li Qiygjj that the formation of this structure depends upon cosmic 
expansion, models of the dark energy in the Universe would be constrained by the 
observation of the large-scale structure. 

One theoretical prediction for the formation of the structure in the Universe, 



. nonlinear perturbation theory, has been considered for a long time. In particular, we 

\ notice Lagrangian perturbation of the cosmological fluid in Newtonian cosmology. 

■ The behavior of the cosmological fluid is constrained by the assumption of a homoge- 

neous and isotropic Universe in the model (e.g., the Priedmann-Lemaitre-Robertson- 
I Walker Universe model). The cosmic expansion is given by Friedmann equations or 

similar equations. Zel'dovich proposed a model of the cosmological fluid based on 
Lagrangian perturbation.® Although the perturbation remains in the linear region, 
density fluctuation in the quasilinear region can be described.lS'lS'lS Once these 
descriptions are available, higher-order perturbative solutions can be derived. Pertur- 
^ . bative equations up to the third order describing dust (pressureless fluid) and a poly- 

tropic fluid have previously been developed^ ^^'^^'^^^'^^'^^ 
An improved method was subsequently proposed j'^^'EB in which the procedure for 
determining the expansion of the perturbation was changed; however, a signiflcant 
improvement was not seen. 
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Recently, resummation theories have been proposedj^''^''^''^''^ these theo- 
ries reorganize an infinite series of diagrams into the form of standard perturbation 
theory and partially resum them. Resummation theories can precisely describe phys- 
ical quantities such as power spectra and correlation functions. In particular, theories 
based on Lagrangian perturbation have several merits: they make is easy to describe 
physical quantities in not only real space but also in redshift spacd^''^''^ and non- 
local bias is included .'^''321' A formula for arbitrary-order (m-loop) correction has 
been proposed.® If we have greater knowledge of higher-order Lagrangian pertur- 
bations, we can further improve resummation theories that are based on Lagrangian 
perturbation. For example, recently, the matter bispectrum with higher-order La- 
grangian perturbations was computed .123 '133 

Although perturbation theory cannot describe strongly nonlinear regions, it still 
plays an important role. Cosmological A^-body simulations have been applied to 
describe the evolution of strongly nonlinear regions. The initial condition for the 
simulations is given by Lagrangian perturbation. It is known that the accuracy of 
the initial condition affects time evolution in strongly nonlinear regions.'3Sl'''3Zli'E3''3Sl' 
In the future, this knowledge would be important for precision cosmology. 

The application of Lagrangian perturbation theory (LPT) to observational re- 
sults will also become important. For example, the Sloan Digital Sky Survey (SDSS) 
Data Release 7 (DR7) contains spectroscopic data on approximately one million 
galaxies.!^ In other words, the SDSS project provides a three-dimensional map 
of the large-scale structure. Because the magnitude limit of the project is approxi- 
mately 22, the survey is carried out within the low-z region. In future surveys, high-z 
{z > 1) galaxies are expected to be observed. It is expected that the results of these 
surveys will provide us with knowledge on the evolution of the large-scale structure. 
Furthermore, because the high-z region is still quasi-nonlinear, higher-order LPT will 
be valid. We expect that very precise predictions regarding the large-scale structure 
can be made by incorporating fourth-order perturbation into resummation theory. 

This paper is organized as follows. In Sec. [21 we briefly introduce the basic 
equations for Lagrangian perturbation; in this study, we consider only a dust fluid. 
Next, we present the perturbative equations in Sec.[3l The fourth-order perturbative 
equations consist of six longitudinal-mode (scalar) equations and four transverse- 
mode (vector) equations. We also derive the perturbative solutions in the Einstein- 
de Sitter (E-dS) Universe model. Every mode grows at the same rate; thus, it 
is not possible to ignore any part of these modes. Next, the application of these 
equations to our results is considered. In Sec. HI we discuss the applications of 
the higher-order perturbations. In Sec. 14.11 we describe one possible outcome of 
applying our result to resummation theories. In Sec. 14.21 we consider initial condition 
problems in cosmological A^-body simulations. Finally, we summarize our flndings 
in Sec. [SJ In Appendix El we describe the relation between the determinant of the 
displacement matrix and the combination of triplet terms. Appendix iBl presents a 
detailed derivation of the fourth-order perturbative equations. Appendix [Cl presents 
higher-order loop correction for the power spectrum. 



4-LPT for cosmological dust fluid 



3 



§2. Basic equations 

In this section, we briefly introduce the Lagrangian description of a Newtonian 
cosmological fluid. When the scale of objects is smaller than that of the horizon, this 
description is valid. The effect of cosmic expansion is reflected by the scale factor a. 
The solution a is given by the Friedmann equations (or alternative equations). In 
the comoving coordinates, the basic equations for a cosmological dust fluid are 

^ + iv,.{^(l + J)} = 0, (2-1) 
ot a 

dv 1 , „ . d 1 ^ , , 

— + -{v-V^)v + -v = -g, 2-2 
ot a a a 

X g = 0, (2-3) 
V^g = -4TTGpba6, (2-4) 

6^^. (2-5) 
Pb 

In Eulerian perturbation theory, the density fluctuation 5 is regarded as a perturba- 
tive quantity. However, in LPT, displacement from a homogeneous distribution is 
considered: 

x = q + s{q,t), (2-6) 

where x and q are the comoving Eulerian coordinates and the Lagrangian coordi- 
nates, respectively, s is the displacement vector, which is regarded as a perturbative 
quantity. From Eq. (|2-6p . we can solve the continuous equation (j2-ip exactly. Then, 
the density fluctuation is given in the formally exact form 

6 = 1- J-\ J = det fl^V (2-7) 



J refers to the Jacobian of the coordinate transformation from Eulerian x to La- 
grangian q. Therefore, when we derive the solutions for s, we can determine the 
evolution of the density fluctuation. 

The peculiar velocity is given as follows: 

V = as. (2-8) 

Then, we introduce the Lagrangian time derivative 

d a 1 , , 

at ot a 



Taking the divergence and rotation of Eq. ()2-2p . we obtain the evolution equations 
for the Lagrangian displacement: 



• ( s + 2-S ) = -47rGp5(J"i - 1), (2-10) 



a 



a 



X ( S + 2-S ) =0, (2-11) 
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where Q refers to the Lagrangian time derivative (Eq. ()2-9p ). To solve the perturba- 
tive equations, we decompose the Lagrangian perturbation into its longitudinal and 
transverse modes: 

s = VS + s^, (2-12) 
V ■ s'^ = 0, (2-13) 

where V refers to the Lagrangian spatial derivative. 
Here, we expand the Jacobian: 

J = 1 + V^S + \{{V^Sf - S,i,S,,, - sl^sl, - 2S,,sl,] 

+ det (% + s^j) . (2-14) 

Eqs. (j2-10p and ()2-lip include the Eulerian spatial derivative, and we convert it to 
the Lagrangian spatial derivative: 
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d 


d 


d 




dqi 


dqj 






d 


d 


d 




dqi 




~\~ Sj iSj^ j 

oqk 



+ ■■■ . 



§3. Lagrangian perturbative equations 



3.1. First-order perturbative equations 

The first-order perturbative equation in the longitudinal mode was derived by 
Zel'dovich.'^'' The evolution equation is given as follows: 

^(1) + 2-5(1) - ATrGpbS^'^ = 0. (3-1) 
a 

The solution can be divided into its spatial and temporal parts: 

5«=5i(t)V'^')(q), (3-2) 

91 + 2-51 - ATrGpbgi = 0. (3-3) 
a 

In the E-dS Universe model (a(t) oc t'^^'^), the solutions to the temporal part are 
given as 

5i(t)oct2/3,t-i. (3.4) 

Hereafter, we consider only the growing solution {gi{t) = t^^^). The spatial part 
■0^1) (q) is given by the initial condition. 

The effect of the transverse mode was examined by BucherlP^ and Barrow and 
SaichPJ 

Vx (■^m+2U'^^A =0. (3-5) 
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We choose a reasonable boundary condition for this equation. The spatial differential 
operator can be removed, leaving 

^ni)+2^^^«=0. (3-6) 

a 

In the E-dS Universe model, the solution is given as 

s^«oct0,t-^/3. (3-7) 

Because there are no growing solutions for first-order perturbation in the transverse 
mode, we consider only the longitudinal mode in the first-order perturbation. 

3.2. Second-order perturbative equations 

The second-order perturbative equation in the longitudinal mode was derived 
by Bouchet et al^^ and Buchert and Ehlers.'^ The evolution equation for the 
second-order longitudinal mode is derived from Eq. (|2-10p . This equation is given as 
follows: 

+ 2^5g) - 4.Gp,S^ = -2.Gp, {s^ S^] - S^]) . (3-8) 

We decompose the perturbation into its spatial and temporal parts: 

5(2) = <^2(0V'('^(q). (3-9) 
The evolution equation is then given as follows: 

92 + 2^92 - ^TrGpm = -AnGpbgl (3-10) 

^S^ = ^KVg-^«^«}. (3-11) 

In the E-dS Universe model, the special solution for the temporal part is given as 
follows: 

<72(t) = -^t'/'. (3-12) 
The equation in the transverse mode is given as follows: 

The evolution equation for the second-order perturbation coincides with that of the 
first-order perturbation. In other words, there are no perturbative solutions for the 
second-order transverse mode. 

3.3. Third-order perturbative equations 

The third-order perturbative solutions in the longitudinal mode were derived 
by Buchert,'i3} Bouchet et al.^^ and Catelan.'^S' The evolution equation for the 
third-order longitudinal mode is derived from Eq. (|2-1U|) : 
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-51 (52 - gl) {V'gVS - ^'S^^l} - ^af det [^pl 



(1) 



(3-14) 



Here, we substitute the first- and second-order equations, ()3-ip - ()3-3p and ()3-9p - 
p-lip . respectively, into the third-order equation. 

The transform between the determinant and the triplet terms is described in 
Appendix |Al We decompose the perturbation into two spatial and two temporal 
parts: 

S^'^ = gsam^^'^Kq) + g3bm^'''Hq). (3-15) 
The evolution equations are given as follows: 



53a + 2-53a " 4:TTGpbg3a 

a 

hb + '^-gsb - ^T^Gpbgsb 
a 



-8TrGpbgi{g2 - gf), 
-MGpbgl, 



(3-16) 
(3-17) 



e^ = ^{^gVS-^g^g}, (3-18) 

= - + (3-19) 

In the E-dS Universe model, the special solutions to the temporal parts are given as 
follows: 

53a(t) = ff, (3-20) 

g?At) = -\t''- (3-21) 

Sasaki and Kasaff^J derived third-order perturbative solutions to both the longi- 
tudinal and transverse modes. In this study, we only consider the effect of the first- 
order longitudinal growing mode. The evolution equation for the third-order trans- 
verse mode is derived from Eq. ()2-lip . By using first- and second-order longitudinal- 
mode equations, (|3-ip - (|3-3p and (|3-9p - (|3-lip . respectively, the equation for the third- 
order transverse mode is given as follows: 

^^,k (slf + 2^<f ) - 51 (52 - g^)e^Jk^^■^^PM " 5i52£^,A:V'SVS = 0- (3-22) 
We decompose the perturbation into its spatial and temporal parts: 



.P=53T(t)Cf(9). 



(3-23) 
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The evolution equations are given as follows: 

gsT + 2^g3T = 47rG/9fe5?, (3-24) 

-^Cf = (^SVS-^>SO,- (3-25) 

Even if we consider only the longitudinal mode in the first-order perturbation, a 
third-order perturbation appears in the transverse mode. The meaning of the trans- 
verse mode will be examined in Sec. [SJ 

In the E-dS Universe model, the solution for the temporal part is given as follows: 

g-Mt) = ^t^. (3-26) 

3.4. Fourth-order perturbative equations 

Although fourth-order perturbative solutions have been derived in the past (for 
example, AijiiiH^^ these solutions were reported before the publication of. -A' There- 
fore, these derivations ignored the effect of transverse modes. In this paper, we 
consider the transverse modes. 

The fourth-order longitudinal mode is also affected by the third-order transverse 
mode, which is generated by the first- and second-order longitudinal modes. 

The derivation of the fourth-order perturbative equations is presented in Ap- 
pendix [Bj The fourth-order perturbative equations for the longitudinal mode can be 
derived. In this section, we decompose the perturbation into twelve different parts. 



six spatial and six temporal: 

+9Mm^'''Hq)+94em^''\q)+9,f{t)4'^\q). (3-27) 
The evolution equations are given as follows: 

<74a + 2-g4a " ^T^Gphg^a = -SvrGpfcgi (gsa - 2gig2 + 25-1), (3-28) 
a 

CL n 

94b + 2-546 - A7rGpbg4b = -8vrGp6gi(c/35 - 2g^), (3-29) 
a 

94c + 2-g4c - ^TrGpbg^c = ^irGpbgf, (3-30) 
a 

9M + 2-gAd - ^nGpbgid = -4vrGpb5(2(52 -^gf), (3-31) 
a 

he + 2-54e - ^TrGpbgAe = ^7rGpbgj{2g2 - gf), (3-32) 
a 

54/ + 2^54/ - AnGpbgif = AnGpbgf, (3-33) 
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^r=^S^Cg\ (3-36) 
^y^ = ^(^SVS-^g^^g), (3-37) 

^i^^ - y^m^^ + ^^i^- (3-39) 

In the E-dS Universe model, the special solutions to the temporal parts are given as 
follows: 



54a (t) = 


(3-40) 




(3-41) 




(3-42) 


gUt) = 


(3-43) 


<74e(t) = -^t'/', 


(3-44) 


54/(t) = 


(3-45) 



Next, we consider the evolution equation in the transverse mode. The evolution 
equation for the fourth-order transverse mode is derived from Eq. ()2-lip . We can 
decompose the perturbation into eight parts, four spatial and four temporal: 

= 54a(t)Cr^(g) +546(t)Cr(9) + <?4c(t)Cr^(g) +54d(t)Cr^(q). (3-46) 

The fourth-order perturbative equations for the transverse mode are decomposed as 
follows: 

QiTa + '^-QATa = -87rGpbgl{g2 - gf), (3-47) 
a 

hrb + 2^54Tf, = -SvrGpfegi , (3-48) 

hTc + 2-giTc = 4^Gpf,5i(53T + gf), (3-49) 
a 

gATd + ^-g^Td = -"^T^Gpbgj, (3-50) 

a 

-vM"' = (*!.;'v'2r'-*2*S''')v (3.51) 
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-V\f' = (*!,;'4?-^.SC.'f)^^, (3-53) 
-V^C.'"' = - 'P^€l€f) ■ (3-54) 

In the E-dS Universe model, the special solutions to the temporal parts are given as 
follows: 

94Ta{t) = (3-55) 
94Tb{t) = (3-56) 
94Tcit) = (3-57) 
94Td{t) = -^i'^'- (3-58) 

The fourth-order perturbative equations in this subsection are the main result 
of this paper. In E-dS Universe model, is proportional to a*^. 

For the generic Friedmann Universe model, we must solve the equations for the 
temporal parts. According to the behavior of 52 and g^, 54 would be proportional to 
mostly. In the matter-dominant era, the solutions almost coincide with those in 
E-dS Universe model. 



§4. Implication of fourth-order solutions for practical applications 

4.1. Applications to resummation theory 

Although the precision of the perturbation improves with increasing order, the ef- 
fect of higher-order perturbation is quite small. When the density fluctuation evolves 
in the nonlinear region, it is difficult to describe the evolution using perturbation the- 
ory. Recently, a set of remarkable methods that substantially improve the precision 
has been studied. These methods constitute the resummation theory.!^ ''^''^^'EU 

There are many methods of the resummation theory, and still new methods have 
been also developed. Hiramatsu and Taruy^^ have developed numerical algorithm 
to solve closure evolutions in both perturbative and nonperturbative regimes. For 
time evolution of matter power spectra, the calculation result in nonperturbative 
regime almost coincides to that from A^-body simulations. Crocce, Scoccimarro and 
Bernardeail^-' have developed fast renormalized perturbative scheme. The feature 
of the scheme is to apply multi-point propagators. They can derive nonlinear power 
spectra quickly. PietronP^ proposed a different method. The method is similar to 
the familiar BBGKY hierarchy.'^ To solve group of evolution equations, we must 
introduce some kind of truncation. In this paper, as one assumption, a part of four- 
point function is truncated. In other words, four-point correlator is consist from 
two-point correlators. Of course, these works consider generic Friedmann Universe 
models. In these works, although the growth rate of n-th order perturbation is given 
by g", the power spectrum is calculated with high precision. 

As we mentioned in Sec. O the relation between the density fluctuation and the 
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Lagrangian perturbation is nonlinear. In other words, Lagrangian linear perturba- 
tion includes a nonlinear effect in the density fluctuation (Eulerian perturbation). 
Therefore, if we choose Lagrangian perturbation as the model, a subset of the infinite 
series of diagrams in Eulerian perturbation would be included.'^ Furthermore, the 
resummation based on LPT can compute the power spectrum in redshift space. The 
mapping from real space and redshift space is given analytically. 

Recently, Okamura et a/.^' proposed a generic formula for m-loop resumma- 
tion via LPT. The formula is shown in Appendix O Their formula can apply to 
generic Friedmann Universe models. They described the evolution of baryon acous- 
tic oscillation (BAO).'^''^ Additionally, they compared improved LPT to A^-body 
simulations in real space. Their two-loop correction can predict the power spectrum 
with a higher accuracy than when using a one-loop correction. However, because 
one-loop correction is sufficiently accurate for describing the nonlinearity, the effect 
of the two-loop correction is shown to be less significant for two-point correlation 
functions. 

Their generic formula is applied to not only the power spectrum in real space but 
also in redshift space. Furthermore, their formula is adopted in the generic Fried- 
mann Universe model. In this paper, we show the apparent form of the perturbative 
solutions in a flat Universe model without a cosmological constant. 

When we consider the transverse mode in third-order LPT (3LPT), at least the 
power spectrum in redshift space would change. The transverse mode in 3LPT does 
not affect the growth of the density fluctuation ()2-7p . However, the displacement 
whose direction is different from that of the longitudinal mode is given. Details of 
the difference between the directions of the longitudinal and transverse modes will 
be discussed in next subsection. 

For application of the resummation theory, higher-order perturbative solutions 
are the point. For two-loop correction for the power spectrum, we require the per- 
turbative solutions up to 3LPT. For more improvement of the power spectrum, we 
will consider three-loop correction. As we show in Appendix [Uj 5LPT is required for 
three-loop correction. 

Making a comparison between one- or two-loop corrections and A^-body simula- 
tions, one finds that two-loop LPT is valid for k < 0.3h~^ Mpc~^ at z = 2 with 1% 
accuracy. 

On the other hand, when we consider bispectrum, 4LPT solutions are enough to 
compute at one-loop order in RPT .EH'-EBxhe contribution of the one-loop correction 
is closed by 4LPT. 

Given the effects of higher-order LPT solutions, RPT offers a more effective 
prediction method for future deep-space surveys. 

4.2. Initial conditions for cosmological N-body simulations 

Although Lagrangian nonlinear perturbations describe quasi-nonlinear evolution 
well, they cannot describe strongly nonlinear evolution. After the formation of caus- 
tics, the hydrodynamical description no longer has physical meaning. For strongly 
nonlinear evolution, cosmological A^-body simulations have been used. With im- 
provements to algorithms for the computation of gravitational interactions and ad- 
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vancements in computer technology, large, detailed simulations can be carried out. 

Although the growth of the density fluctuations starts in the recombination era 
(z ~ 1000), A^-body simulations are begun from a later time, because of the difficulty 
in performing the simulation with small fluctuations. In the standard scheme, the 
small fluctuations evolve following the perturbative approach until they reach the 
quasi- nonlinear regime {z ~ 30) 

In the perturbative approach, ZA has long been used. In ZA, although the 
perturbation remains in the linear regime, quasi-nonlinear density fluctuations can 
be described. However, ZA does not reproduce the higher-order statistics, because 
the acceleration is always parallel to the peculiar velocity. In the quasilinear regime, 
effects that cannot be expressed using only linear perturbations appear. If we set 
up the initial condition for use with a cosmological A'^-body simulation at the quasi- 
nonlinear stage, ZA seems insufficient for describing the initial condition. 

To solve this problem, Crocce, Pueblas, and Scoccimarro^''' proposed an im- 
provement in which different initial conditions are adopted. Essentially, the ini- 
tial conditions are based on approximations that are valid up to second-order LPT 
(2LPT) and reproduce the exact value of the skewness in the weakly nonlinear re- 
gion. They set up the initial conditions for ZA and 2LPT at z = 49. The simulations 
were run using the Gadget2 cod^® and contained N = 512^ particles. The cosmo- 
logical parameters used were l?m = 0.27, f2y\ = 0.73, f2b = 0.046, h = 0.72, and 
asiz = 0) = 0.9. The box size was Lbox = 512, W2Ah~^ [Mpc]. 

Crocce, Pueblas, and Scoccimarro showed the evolution of higher-order cumu- 
lants such as skewness and kurtosis. First, they computed these quantities for the 
reference runs (2LPT, Zi = 49). Next, they compared the quantities for reference 
runs with those for several simulations performed using other initial conditions. The 
difference between the skewness and the kurtosis values for the reference runs and 
the simulations using ZA is several percent for 2 < R < 60h~^ [Mpc] at ^; = 0, 1, 2, 3. 
Simulations using ZA underestimated the power spectrum in the nonlinear regime 
by approximately 2%, 4%, and 8% at z = 0, 1, and 3, respectively. 

Thereafter, the effect of 3LPT was verified by Tatekawa and Mizuno. The condi- 
tions they considered differed from those in the former study. The simulations were 
run using the P^M code^''' and contained N = 128'^ particles. The cosmological 
parameters were taken from the WMAP three-year results.'^ They considered only 
the longitudinal mode in 3LPT. Their initial condition was given by z ~ 23, 33, 84. In 
their results, the difference of both the skewness and the kurtosis appeared between 
the simulations using ZA and 2LPT. Even if the initial condition is given at z ~ 84, 
the difference appeared. The result is consistent with that by Crocce, Pueblas, and 
Scoccimarro. On the other hand, if the initial condition was given by z ~ 33, the dif- 
ference between both the skewness and the kurtosis for the simulations using 2LPT 
and 3LPT is less than 2% for R = 2h^^ [Mpc] at late times. According to their 
results, 2LPT is sufficiently accurate when several percentage points of precision are 
needed for the higher-order cumulants. In other words, the initial conditions given 
by 2LPT have the enough accuracy in simulations with precision of one percent. 

In the future, when higher precision is needed, we should select the initial con- 
ditions for A^-body simulations using 3LPT or 4LPT. The transverse mode of 3LPT 
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was ignored in previous analyses. Before estimating the effect of 4LPT, it is necessary 
to estimate the effect of the transverse mode in 3LPT. 
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Fig. 1. The displacement generated by Eq. H4-l|) . In this case, ip^f^'' disappears. Here we plot in 
the range of §1,52 £ [0, 7r/4]. In these figures, the vectors are amplified, (a) ILPT; (b) 2LPT; 
(c) 3LPT longitudinal mode; (d) 3LPT transverse mode. 



Here we suppose simple perturbation in first-order LPT (ILPT): 

7/^1(9) = e {cos gi + cos ^2} • (4-1) 

This perturbation gives the concentration of the matter. The displacement vector is 
shown in Fig. [TJ The longitudinal mode in 2LPT and 3LPT promotes the concen- 
tration. The direction of the displacement in both 2LPT and the longitudinal mode 
in 3LPT is parallel to that in ILPT. However, the displacement of the transverse 
mode in 3LPT differs from that of the others. This difference in the direction of the 
displacement would especially affect the statistical quantities in redshift space. 

Here, we roughly estimate the effect of 4LPT. The difference between the higher- 
order cumulants in the cases of ZA and 2LPT was approximately 10%. Li contrast, 
the difference between the higher-order cumulants in the cases of 2LPT and 3LPT 
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was approximately 1%. When the order of the perturbation is raised, the difference 
is 1/10. This tendency suggests that when the desired precision is less than 1%, we 
should select the initial condition using 4LPT. 

§5. Summary 

We have derived Lagrangian fourth-order perturbative equations for a cosmo- 
logical fluid. The equations consist of six scalar equations and four vector equations. 
The spatial components of the perturbation are described by Laplace's equation. 

In this paper, although we assumed that the linear perturbation only had a 
longitudinal mode, a transverse mode appears in the third-order perturbation; we 
consider the effect of this transverse mode. The density fluctuation is described 
by Eq. (j2-7p . In 3LPT, even if the transverse mode is present in the third-order 
perturbation, the mode does not affect the density fluctuation, because the mode 
disappears in Eq. (|2T4p : this mode affects only the displacement. The simple case 
was shown in Sec. 14.21 

In 4LPT, the doublet term consists of ZA, and the third-order transverse mode 
appears. When the order of the approximation is raised to fourth order, it is neces- 
sary to consider the effect of the transverse mode on not only the displacement but 
also the density fluctuation. 

The appearance of the transverse mode does not imply the occurrence of vortic- 
ity. In LPT, the vorticity is described exactlyPJ'IS'lSJ 

u}{q,tQ) refers to the initial vorticity. Therefore, if vorticity does not exist in the 
initial condition, it never appears. This result is consistent with Kelvin's circulation 
theorem. 

Here, we note a problem in higher-order perturbation. If the primordial density 
fluctuation is negative, the fluctuation 6 approaches —1. In Eulerian perturbation, 
the density fluctuation exceeds —1, i.e., negative mass density appears. However, in 
LPT, 6 decreases gently and converges to —1. However, the first-order perturbation 
remains the best approximation for describing the late-time evolution of voids. In 
particular, if we stop expansion when an even order (second, fourth, sixth, . . . ) 
is reached, the perturbative solution describes the contraction of a void.'^''^''^ 
Therefore, when significantly advancing the time evolution, we should apply the 
fourth-order perturbation. 

Resummation theory significantly improves the precision of the perturbation. 
Even if we consider fourth-order Lagrangian perturbation, we could still avoid the 
above problem. In particular, the Lagrangian perturbation easily computes physical 
quantities not only in real space but also in redshift space. In the future, the resum- 
mation of the Lagrangian perturbation up to the higher-order will contribute to the 
prediction methods used in precision cosmology. 

Rampf and BucherlPS have been derived 4LPT equations independently. Al- 
though their peculiar-system approach corresponds to our approach, their 4LPT 



14 



T. Tatekawa 



equations are different from our equations. The origin of the difference are treat- 
ment of temporal components and spatial derivative. For temporal components, they 
introduced conformal time for the equations. Then they considered both growing 
and decaying modes. For the spatial derivative, they considered inverse matrix for 
the Jacobian matrix between Eulerian coordinate and Lagrangian coordinate. At 
present, especially in the spatial components, both equations are inconsistent. The 
meaning of the difference in the spatial terms and its effects should be confirmed in 
near future. 

For more higher-order correction, 4LPT would not be enough. Enough higher- 
order correction to resummation, fifth-order LPT (5LPT) has been achieved. Based 
on the calculation of 4LPT, 5LPT equations will likely consist of at least 19 scalar and 
18 vector equations with our approach. Although the derivation of 5LPT equations 
appears to be quite a challenging problem, useful predictions for future deep-space 
observation might be possible; therefore, it is worth wrestling with this problem. 
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Appendix A 

Determinant of displacement 

We consider rewriting the determinant of the Lagrangian displacement matrix 
with triplet terms. 

First, we define V as follows: 

V {A, B, C) = 22C,33 + Ai2 -8,23(^,31 + ^,13-6,21(^,32 

—^,11-6,23(^,32 + ^,22-B,3lC 13 -I- ^,33-8,12(^,21 
+^,33-B,llC,22 + ^,31-6,12(^,23 + ^,32-B,l3C,21 
— ^,32-B,ll(7,23 + ^,13-6,22(^,31 + ^,21-6,33(^,12 
+^,22-6,33(7,11 -I- ^,23-6,31(7,12 + ^,21-6,32(7,13 
— ^,23-6,32(7,11 -I- ^,31-6,13(7,22 + ^,12-6,21(7,33. 

= ^,n-6,jj(7,fcfc -|- A^ijBjif;C^ki + A^ijB^iiiCjk 

V {A, B,C)=V {B, C,A)=V (C, A, B) 
= V {A, C,B) = V {B, A,C)=V (C, B, A) , (A-2) 

V{A,A,A) = 6 det {A , (A-3) 

where A, B, and C are scalar functions of q. In this paper, we consider fourth-order 
perturbation. Therefore, we must consider a determinant that consists of first- and 
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second-order perturbations: 

det (sj) + = (5(1) + 5(2), 5(1) + 5(2), + 5(2) 

= ^[v (5(1), 5(1), 5(1)) + p (5(2), 

+0(e') 

= (5(1) , 5(1) , 5(1)) + (5(1) , 5(1) , 5(2)) + 0(e5). (A.4) 
The determinant is rewritten in triplet terms as 



det 5 



(5(1), 5(1), 5(1)) 



6 

1 0(1) 0(1) 0(1) 1 qW c(l) c(l) , 1 c(l) 0(1) 0(1) /- A c:^ 



(5(1), 5(1), 5(2) 



(1) ^(1) _ ^(.J^; ^ 



(l)c(l)\5(2) 



(1)0(1)0(2)' 



_l_9 / 0(1) 0(1) 0(2) _ qyi-l Q\'-) 
\ Jfc "^,11 '-'jk'^Jk 



(A-6) 



Appendix B 

Derivation of fourth- order perturbative equations 



In this section, we provide a detailed description of the derivation of fourth-order 
perturbative equations. First, we present a derivation of the longitudinal mode. The 
left-hand-side terms of the evolution equation ()2-10p for the fourth-order longitudinal 
mode are given as follows: 



5^^+2-5 



(4) 

,ii 



(3)/5(^)+2-5(^) 



a 



o(l) f .(3) , (3) \ (1) (1) / x(2) a ■ (2) 



, 0(1) 0(2) / ^(1) qW \ , c^(2) o(l) / o(l) , 9« 6(1) 

,ij ,jk \'^,ki "T'^^'^.fci ) '^'^,ij'^,jk y^,ki 



0(1) 0(1) 0(1) 
^,jk^,kl 



(B.l) 



The right-hand-side terms of the evolution equation ()2-10p for the fourth-order lon- 
gitudinal mode are given as follows: 



0(4) ,90(1)0(3) , o(2)c;,(2) 



1/0(2)0(2) 0(2) 0(2) \ 

ro(l)o(3)_ o(l)J3)1_ 1^/^(1) ^(1) ,.(2) 
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(1)5(1) 



o(l)c.(l)\ / 0(1)0(1) 



0(1)0(2) 
jj ^,kk 



+ det ( 5 



^(1) 



o(l) 0(2) 



00(1) 0(1) 0(2) 



0(1) 0(1)" 

0(1) 0(1) 
jj ^,kk 



0(1) q(i) / 0(1) 0(1) 0(1) 0(1) 1 , 0(1) 0(1) 0(1) 0(1) 
2 Jj X^^M ,11 '^,kl ,kl j •^jj'^M M 



0(1) 0(1) 

^,jk^,jk 



(B-2) 



The definition of V is given in Appendix 1X1 We rewrite the homogeneous terms in the 
fourth-order perturbative equation using first-, second-, and third-order equations 
(iMj-dMI), (IMD-dMiD, (I3T6]1 - (I3491 ). and ([3^ . respectively. First, we examine 
the doublet terms in the first- and third-order perturbations: 



(3) / o(l) 



.a ■ 



9q(1)o(3) 



^(1) 



{ 



0(1) q(3) 



:(i) J3) 



(3)1 



3^„/,(l)„/,(3b) , ^4 ,,(1)^(3) 



/ ,(1) , (3a) ,{1) , (3a)\ /"w/^) .A^^'') .A^^^.A^^''^^ 



gig?,. 



-gigsTTp^ljCij 



-87rG/>69i(53a - 25152 + 2^?) • \ (V'SV,' 



(l)„/,(3<^) _,/,(l)„/,(3«) 



-STrGpbgi {gsb - 25' 



',(l)./,(3'') 



',(l)./.(3^') 



+ 47rGpbgtij^^CS 



+4.Gp,gl{-2g, + 2,?) • \^ {^g^g - ^«^(^} 



-SvrGpbstV'S^det (v,^) 



(B-3) 



The upper two lines on the right-hand side of Eq. (|B-3p correspond to the source 
terms of Eqs. (|3-28|) - (|3-3U|) and Eqs. (|3-34p - (|3-36p . The other hues (remainder terms) 
in Eq. (|B-3p will be considered later. Next, we examine the doublet terms of the 
second-order perturbations: 



S 



(2) 



^TlGpb 



0(2) o(2) 
,ii ,jj 



0(2) r,(2) 
,ii ,jj 



0(2) 0' 
,ij 



(2)1 
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The first line on tlie riglit-hand side of Eq. ()B-4p corresponds to the source terms of 
Eqs. (j3-3ip and p-37p . The last terms (remainder terms) in Eq. ()B-4p will be consid- 
ered later. We notice the triplet terms in the first- and second-order perturbations: 



{2)c(l) / dW , n^Mi) 



,ij ,jk \ ,ki 



.ki 



,ki 



_C.(2) (q{l)q{l) _ q{l) q{l)\ ,00(1)0(1)0(2)" 

,ii \ ,jj ,kk '-'jk'-'jkj ^ '^,jj'-',kk 



AirGph X 



'.(2)„,,(1)./,(1) 



',(2).,.(1)„,,(1) 



(B-5) 



Here, we add the remainder terms in Eqs. ()B-3p and ()B-4p . The triplet terms become 

-i^2^SH^S^S-^S^!^} 

-(2.2-.?)^«^^^g+.2^ff^S^g 



47rGpb5?(252-5?) 



(B-6) 



The first line on the right-hand side of Eq. (|B-6p corresponds to the source terms 
of Eqs. (|3-32p and (|3-38p . The last terms (remainder terms) in Eq. ()B-6p will be 
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considered. Finally, we consider the quartet terms in the first-order perturbations: 



0(1) 0(1) 0(1) 



,(,6 ,1,6 



+4TTGpb X 



gil) ^(1) 



} {s[,:ls 



(1) 

kk'-'Al 



} 



25« det f5« 



,3 (1) (1) ( Mi)qW 

"•"2 .ii \'-^,A:fe'-^,i« 



c.(l)o{l)\ , 0(1)0(1)0(1) 0(1) 
'-',fc« •-',kl J "T "-".ij '^,kk'-',ll 



(87) 



Here, we add the remainder terms of Eq. (|B-6p . The quartet terms are rewritten as 
follows: 

+^*,i'*g{'/'S'/'!.;'-*S'/'S}" 



47rGpbgt x 



(B.8) 



The right-hand side of Eq. ()B-8p corresponds to the source terms of Eqs. ()3-33p and 
(13^ . 

Next, we derive the equations that describe the transverse mode. The evolution 
equation is derived from Eq. ()2Tip as 



^ijk [ ^kj + a^kj '* 
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0. (B-c 



Three terms on the left-hand side of Eq. ()B-9P disappear owing to the symmetry of 
subscripts in the longitudinal mode, and so 

^^^kSfi [s^^ + 2^^g) = 4vrG/>,9i%;^V',?V;3 = 0, (B-IO) 

^^^kSf,S% {S% + 2^5« ) = 4vrGp,5?52e.,.V',^,Vf;i^?j. = 0, (B-ll) 

^^^kS^Ss%S% [S% + 2^5«) = 4vrGp,,^.,.^,(;)^;;i^«^,« = 0.(B.12) 
Equation (|B-9P is rewritten as follows: 

^^'fej ^ + 2~*fc,j ~ ^ijk^Jl (^^kj + ^Q^i,; ) ~ ^ijk^^jj (^^^kl + 

By using first-, second-, and third-order equations ()3-ip - (|3-3p . ()3-9p - (|3-lip . (|3-16p - 
(|3-19p . respectively, and (|3-23p . the equation is rewritten as 



a 



p. ., s ^ ' -\- 2-s ^ 



= ATTGpbSijk X [51 {93b - 251(52 - gf)} ^fi^^^H + 9i{93b - 2gl)^pf,^f^^ 

4 , (1)^(3) I „ „ 

+ amri^^lkt! - . (B-14) 

Upon decomposition of the fourth-order perturbation by Eq. p-46p . we obtain Eqs. 
(CTD - (I3^ . 

Appendix C 

Higher- order loop correction for power spectrum 



The generic formula about the power spectrum with A^-loop corrections in both 
real space and redshift space has been derivedPS Here we suppose Gaussian initial 
fluctuation. In real space, the power spectrum is described as 



P{k) = exp 



2 hi • ■ ■ ki2n A (2n) 

^ (2n)! ^1-^2" 



71=1 



00 



m=l 
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where 



(2n) 



00 / 00 
n=l 



(2n)! 



(C-l) 



■ 2n 

n 



1^ (&^) ,P2n). (C-2) 



(5/5 is Dirac's delta function. Ci^...j2„ is a 2n— order cumulant of the displacement 
field in Fourier space iP: 



•2n, 



,P2r. 



(C-3) 



For derivation of the power spectrum, A^-loop contribution to the power spectrum 
predicted by SPT is required. In A^-loop correction of Lagrangian resumption, the se- 
ries in the exponent is expanded up to O [{Pl)^) ■ The remaining factor is expanded 
up toO((PL)^+i). 

Following the formula (jC-ip . three-loop correction for the power spectrum is 
derived. Here we extend correction terms. 



,(2) 



,(4) 
ijlm 



ij ' ij ' ij ' ij ' ^] 



I (22) 



,(13) 



I (33) 



(24) 



+ 2A. 



(15) 



H 



6^' 



(1122) 



(C-4) 

\jlm ~'~ "^"^ijlm ' ^) 

The superscripts of A^°'^'^ mean that A consists of a-order and /3-order perturbative 
quantities. In Gaussian initial fluctuation, ^(^^^m) does not appear. The explicit 
form of three-loop correction for the power spectrum is shown as follows: 

--^ ( ^(11) + >l(22) + 2^(13) ^ ^(33) ^ 2^(24) ^ 2^(15) 

kikjkikm /3_^(ii22) ^ 2^(^^^3) 

-loi 
SPT 



P{k) = exp 



g \^~" 'ijlm ' ijlm 

PUk) + P^p^r^lA:) + P^pl°°P(fc) + P^p^°°P(^) 



SPT 

U2 



+ p^i^k) I (^(") + ^(2^) + 2^(13) ^ ^(33) ^ 2^(24) + 2^(1^) 



+ 
+ 



k^ 



6~ 



+ 2^(^1)^(22) ^ 2^(11)^(13)) + 



k^ 



12967r6 



(-4( 



11)n3 



o /l(1122) I 2 /I, , . 
"^•^ijlm ^ ■^•^ijlm I 



(1113)^^ 



ik){A' 



P^-l°°P(A:)^(ii) 



67r2'SPT 



(C-6) 



For the computation for ^(i^), 5LPT solutions are required. Therefore, three-loop 
correction for the power spectrum, 4LPT is not enough. 
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